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The phenomenology of Josephson tunnel junctions between unconventional superconductors is 
developed further. In contrast to s-wave superconductors, for d-wave superconductors the direction 
■ dependence of the tunnel matrix elements that describe the barrier is relevant. We find the full I-V 

0^ ' characteristics and comment on the thermodynamical properties of these junctions. They depend 

sensitively on the relative orientation of the superconductors. The I-V characteristics differ from 
the normal s-wave RSJ-like behavior. 
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The symmetry of the superconducting order parameter of the high-T^ compounds is the subject of a heated 
I debate. Quite early there were theoretical suggestions that these materials are unconventional superconductors 
that have da;2_y2-symmetry of the order parameter. They have stimulated both experimental and theoretical work. 
The experimental situation is still unclear 15 1. Several experiments indicate d-wave pairing, e.g., the temperature 
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, dependence of the penetration depth A(T) |] and the NMR relaxation rates The most convincing set of 



JID loop consisting of a high-T, 
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' experiments up to now are the SQUID experiments |9|jlO[|, which show that a SQ 
, crystal closed by an ordinary s-wave superconductor shows a phase shift of the order of tt in the dependence of 
the critical current on the flux, and the experiments on flux quantization in multi-domain rings ||ll|,|l2|. There are 
[ other experiments, however, that cannot be reconciled with the d2,2_j^2 -model and have been interpreted to indicate 
. (possibly anisotropic) s-wave pairing |l^,|l^, see, however, Ref. |l5). Since a microscopic theory is still lacking, there 
' are also different theoretical opinions on the symmetry of the order parameter. Some phenomenological calculations, 
, assuming d-wave pairing, have been done and they yield partial understanding of some experimental data. 

a In this note we will further develop the phenomenology of unconventional superconductors, assuming BCS-like 
I I behavior for the density of states. Our results are not specific for high Tc materials, but may be relevant also for 
• other unconventional superconductors such as heavy fcrmion compounds. We will consider tunnel junctions between 
^ , two pieces of bulk superconductor as shown in Fig. ^ The superconductor is assumed to have a fc-dependent order 
O ■ parameter of the form = Ao[fc^ — ky]. In contrast to earlier work |p^-p^, we include the direction dependence 
' of the tunnel matrix elements that describe tunneling across the insulating barrier between the superconductors. For 
^ . s-wave superconductors this dependence drops out, but for d-wave superconductors it is essential pO| , ^l| ]. 

Along the lines of Ref. j22| we calculate the full current-voltage characteristic (I-V curves) for different relative 
orientations. We find that the presence of gap-nodes strongly influences the I-V characteristics. For a specific relative 
orientation the quasi-particle current will be proportional to the square of the voltage, which leads to new behavior 
for the total current, which is different from the usual RSJ-like overdamped junction. We will also comment on the 
thermodynamic properties ||23|| of d-wave tunnel junctions and the relevance of our work for thin granular high-Tc 
films 

A convenient starting point is the tunneling Hamiltonian 

H^HL+Hn+[ I ^[T(r,r')<<,Vr'.. + M- (1) 

JreL Jr'eR „ 

Here Hl and Hji denote the unperturbed Hamiltonians of the left and right superconductor. From a second order 
perturbation expansion in the tunnel matrix elements one finds the Ambegaokar-Baratoff formula that expresses 
the current across a tunnel junction in second order perturbation theory in the tunnel matrix elements as 
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Here A and i? are the spectral densities of the normal and anomalous Greens functions, </? is the phase difference 
between the two superconductors, and Jl/r the Fermi distributions at the chemical potential ^ll/r- The left and 
right chemical potential differ by the applied voltage, — ^r =^ eV. The Tk^k' are the matrix elements that transfer 
electrons from a state k in one superconductor to a state k' in the other. The first term in (||) describes the quasi- 
particle current Iqp and the second term the supercurrent which is proportional to the critical current Icr- The 
spectral densities have the form 



A{k,uj) = vr 



l + ^^)si.-Em+[l-^^]Si. + Eik)) 



B{k, u) = [6{uj + E{k)) - 5{u: - E{k))] , 



where Eik) = y | A(fc) ^pj^g strategy is to take seriously the orientation dependence [g^^. This is necessary, 

since with the standard assumption that the tunnel matrix elements T^ jf are independent of momenta, Eq. (||) yields 
a critical current Icr = after an angular average with = Ao[fc^ — k^]. 

Thus, in order to obtain physical results, the direction dependence of tunneling has to be taken into account. We 
will consider the case of a smooth barrier for which momentum parallel to the barrier is conserved during tunneling. 
A reasonable assumption seems to be to take the tunnel matrix elements T(r, r') non-vanishing only when r and r' 
are both close to the barrier. We take 

T{r,r')=f5^^\r-r')5^^\r,-xo), (3) 

where is the location of the barrier, as indicated in Fig. |l|. This ensures that the momentum parallel to the tunnel 
barrier is conserved, since the Fourier transform of the matrix element is 

Tu,-k' = (2^)'f <5(2)(fc|| _ |)e(fc, . fc^) (4) 

The last factor is a restriction on the direction of tunneling. Before tunneling the electron should move towards the 
barrier, after tunneling away from the barrier. The function /(| kx |) is a weight function that makes tunneling 
perpendicular to the barrier more probable than tunneling parallel to the barrier. If one models the barrier by a 
slab of finite thickness and takes the tunneling probability to be exponentially small in the traversed distance in the 
barrier, one deduces /(| kx |) '~ exp(— 1/ | kx |). Numerical constants may depend on the exact choice for /, but 
principal behavior should not depend on it. In the following we take it to be /(| fc^ |) =| fcx | which corresponds to 
vanishing thickness of the barrier. 

With this explicit form of the tunnel matrix elements, both the quasi-particle current Iqp as well as the critical 
current Icr following from Eq. ^ can be evaluated by numerical integration for different temperatures. For nor- 
malization we calculate the normal state current for Ao=0 with the help of (|^). This allows us to identify the normal 
state conductance Gjq of the barrier as Gjv = 87r^T'^iV(0)^A:^^^/i?x, where A denotes the area of the junction and 
Rk — h/e^ the Klitzing quantum of resistance. The expression for the critical current that we find from Eq. (H) with 

/i; = /ir is 

Icr = 2e / dujduj I —- ; — Nl(lu)Nr{lu ) \ T{Q,il ) | , (5) 

/ J An An u! — U! oj ui' 



where Nl{uj) = N{0) \ lj \ Q{\ uj \ — Ai(il))/ yju"^ — A'^(ri) denotes the (angle-resolved) density of states of the left 
side, Nr{lu') is defined similarly for the right side, dfi = d(j)ddsm6, and Ai/^(ri) — Aq cos[2((/) — (Pl/r)]- The angles 
(f)L and 4>fi determine the relative orientation of the two superconductors. A natural choice which applies to many 
junctions is (j)L = 4>r = 0. The results are shown in Figs. p| and ^) for Iqp and Icr respectively. The magnitude 
and sign of the critical current show a strong orientation dependence, see Fig. ^). 

The quasi-particle current behaves mostly like for an s-wave superconductor Iqp ~ exp{— Agj^/fc^T}, where the 
leading behavior is determined by tunneling from a gap node in one superconductor into the effective gap Ae/ / in the 
other. However, for those special relative orientations for which the gap nodes in the left and right superconductor are 
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parallel, the behavior is different. The quasi-particle current Iqp ~ for voltages 2Ao > eV > fc^T and Iqp ~ V 
for voltages eV < ksT. This can be understood as follows. The dominant contribution to Iqp arises from 'node to 
node' tunneling. The available phase space around the gap-nodes scales with the bigger of eV and T. With the usual 
factor V from the difference in — fp this yields the quoted behavior. 

The total current / may be obtained by integration over time of Kirchhoff 's equation in the phase representation, 
/ — IcR,sai{(j)) + lQp[2eV — h(f>). In this way the average time-derivative of the phase is found as a function of the 
total current through the junction. The result is shown in Fig. ^ for low temperatures kpT ^ Aq and (pL — (pR — 0; 
i.e. the geometry in which node to node tunneling appears and Iqp V'^. The I-V characteristics deviate clearly 
from the well known RSJ behavior. For higher temperatures T of the order of the gap Aq we have Iqp ^ V (see Fig. 
H) and the RSJ-behavior for the I-V curve is recovered. 

We now continue with the discussion of the thermodynamic properties that are described by an effective action for 
the phase difference across the junction. This is especially relevant for small mesoscopic tunnel junctions with a low 
capacitance C [ p5| . From a second order perturbation expansion in the tunnel matrix elements one finds the following 
effective action for the phase difference (p across the junction 

S[ip\=J dr— j +J drdr cos{ )a{T ~ t ) + J drdr cos{ )/3(r - r ) . (6) 

The kernels a and /3 describe quasi-particle tunneling and the Joscphson coupling respectively. They are given by 

^3^, r Gik,T)Gik',^T) 




F{k,T)F{k\ 



Tk,k' r , (7) 



-T 



where the normal and anomalous Greens functions G and F were introduced. The usual approximation, in which 
the Fourier transform of the tunnel matrix element T^^k' is taken to be independent of momenta, leads for s-wave 
superconductors to the standard expressions in the literature, i.e. both a and /3 decay exponentially on the time-scale 
A~^. The same approximation for a d-wave junction will lead to unphysical results. For instance the beta-term is 
found to be zero, as the average over orientations of A^ is zero. 

To remedy this shortcoming, one again has to retain the direction dependence of the tunneling matrix elements (jj). 
The gap nodes have a pronounced effect on the long time behavior of a and (3. If (/i^ = 0^ = 0, the gap nodes in the 
two superconductors are in the same direction (the "node to node" geometry) and the asymptotic behavior of a{T) is 
T^"^. This corresponds to a low frequency behavior ^ uj'^lnuj, which is 'super-ohmic'. An investigation of the oj'^lnoj 
dissipation remains a subject for further study p6t , and may be relevant for the phase diagram and superconductor- 
insulator transition in thin granular high-Tc films p^ . For relative orientations without node-to-node tunneling we 
find exponential decay, as was also found for s-wave tunnel junctions. 

In conclusion we have calculated the full I-V characteristic for a tunnel junction between two unconventional 
superconductors by making a new Ansatz for the direction dependence of the tunnel matrix elements. The temperature 
dependence and a strong dependence on relative orientation is found. For specific orientations the I-V characteristics 
differ from the usual RSJ-like behavior. All of our predictions have experimental consequences and should be verifiable 
using thin-film junctions on bicrystal substrates. 
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FIG. 1. The system we consider consists of two coupled d-wave superconductors. Their orientation is characterized by the 
angles 0z, and 4>r. 
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FIG. 2. The quasi-particle current Iqp as a function of the applied voltage for temperatures T/Ao=l (upper curve), and 
0.1 (lower curve) for the case 4>l = 4>R=0- Here and in the following figures, 7jv = lN{eV = 2Ao). 
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FIG. 3. a) The critical current Icr as a function of temperature for the case 4>l = 4>r=0 (full line). For comparison the 
s-wave result is also shown (dashed line), b) The critical current Icr as a function of the relative orientation (j)R {4>l=0 (full 
line), and 4>l = j (dashed line)) at T=0. 
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FIG. 4. The complete I-V characteristic for a d-wave tunnel junction with (f>L 
also the RSJ result is shown (dashed line). 
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